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MOTION OF A CYLINDER UNDER THE SURFACE OF A HEAVY FLUID* 
; By L. N. Sretensky 


1. INTRODUCTION 

The present work on the theory of the motion of a solid body in a 
fluid having a free surface consists of two parts . 

In the first part (sections 2 to 8), general equations are given 
for the determination of the flow of a heavy fluid of infinite depth 
about a submerged circular cylinder. The problem of the motion of a 
cylinder under the surface of a heavy fluid presents considerable dif- 
ficulties in its solution. These difficulties were first pointed out 
by Kelvin. A solution is given herein for the simplest part of the 
problem of Kelvin, namely, setting up the equations of the problem and 
obtaining certain approximations of their solution. The approximate 
solution obtained replaces the moving circular cylinder by a certain 
vortex . 

T. H. Havelock (reference l) in a recent paper considers the prob- 
lem of Kelvin under the same general assumptions as are herein con- 
sidered, but gives a more advanced approximate solution. 

In the second part of the paper (sections 9 and 10), the cylinder 
is replaced by a dipole of a certain strength, and an equation is set 
up for the computation of the wave resistance of a circular cylinder 
moving in a fluid of finite depth. 


2. DERIVATION OF BOUNDARY CONDITIONS OF PROBLEM 

In order to study the problem of the motion of a cylinder under the 
surface of a fluid, a system of Cartesian coordinates XOY is intro- 
duced. The OX-axis is placed along the undisturbed surface of the 
liquid with its positive direction coinciding with the velocity c of 


"Dvizhenie tsilindra pod poverkhnostyu tyazheloi zhidkosti." 
CAHI Rep. No. 346, 1938, pp. 3-27. 
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the flow at infinity, and the OY-axis is taken vertically upward. The 
radius of the cylinder is denoted hy a; and the ordinate of its center 
by h. The origin of the XOY system is chosen directly above the 
center of the cross section of the cylinder. 


The fluid moves with the velocity potential $ and the stream 
function i|r . The disturbance (due to the presence of the cylinder) of 
the horizontal projection of the fundamental flow velocity c is 
denoted by uj the disturbance of the vertical projection, by v. If 
c is the horizontal velocity of the flow at infinity, the following 
expression can be obtained for the projection of the total velocity v 
of the particles of the fluid: 

a + c, v (!) 

The equation of Bernoulli can now be written, if the motion is assumed 
to be steady: 

T-c-ev-^v' (2) 

where & is the mass density of the fluid. 


Along the free surface the pressure has the constant value Pq. 
Equation (2) is applied to the particles of the free surface lying far 
ahead of the cylinder, that is, to those for which x = - •». For 
x = - °°, V = c and y = Oj therefore. 


Equation (2) can then be written 




c') 


But 


V*' = (« + «}* -{-$* 


hence 


A _ gy - ^ + + 2cu) ( 3 ) 

The increments u and v are now assumed to be so small that their 
squares may be neglected. With this assumption, the accurate equa- 
tion (3) becpmes the approximate equation 

P-Po 


o 


= —gy~cu 
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When this equation is applied to the free surface. 



( 4 ) 


This equation can he used to determine the shape of the disturbed sur- 
face of the fluid. 


The fact that the velocity of the particles moving along the surface 
is directed along the tangent to the surface is now taken into 
consideration: 

dy _ v 
dx~ u-\-c 

This relation may be represented in the form 


d l— v Ti 
dx c [ 



If the terms of the second and higher order smallness are rejected. 


dy v 

dx~ c ( 5 ) 


Eliminating the ordinate y 


by combining equations (4) and (5) yields 


dii 

dx 



( 6 ) 


By introducing the potential <p(x,y) and the stream function \|r(x,y) 
of the disturbed velocities, this condition may be given another form: 


tt. 


d<p 

dx~~ dy 


v = = ( 7 ) 

dy dx 

The functions <p and \|r are harmonic conjugate functions. Substi- 
tuting expressions {l) in condition (6) yields 


' g d’j » 

dxdy ~c* <fx 


( 8 ) 


A second simplification of the problem is now introduced by requiring 
that y be replaced by zero. Equation (8) will then be satisfied 
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along the’ axis of abscissas. Integrating along this axis gives 


( 9 ) 

where the constant of integration is zero from the consideration that 
the zero value of the stream function Y was ascribed to the surface, 
and moreover, for x = - <», u = - cty/dy = 0. 

Condition (9) is the first, boundary condition of the problem. The 
conditions on the surface of the cylinder will now be discussed. Let 
the cylinder be washed by the streamline f = a. Since the relation 
between f and \|r is given by 


»J- = cy-\- ■'< 


the condition Y = a on the cylinder will take the form 


= 1 4 - cy 

where a > 0 . 

The determination of the motion of the fluid thus depends on the 
integration of the Laplace equation 


d'^ 

dx? 



for the following boundary conditions: 


for V — 0 
dy c* ' ' 

•i = a 4 -cy for x- -f (v + hf = o “ 


(I) 

(II) 


3. TRANSFORMATION OF FLOW REGION INTO CIRCULAR RING 

In order to obtain the integral of the Laplace equation correspond- 
ing to conditions (i) and (il), a conformal mapping of the complex var- 
iable plane z = x + iz on the plane £ = £ + if) is carried out by 
setting 
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z — )i 
z \^li 


where 


/. = — V h- — a-< 0 


(i) 


( 2 ) 


When the complex variable z is varied along the OX-axis from -® to 
+« the point £ describes, starting from the point £ = 1, the circle 

£ 2 + |j2 _ 2 _ i n the clockwise direction. 

If z = - hi + ae^ is set next in formula (l) .and the angle y 
is varied from zero to 2«, the surface of the cylinder will correspond 

to the circle £ 2 + i} 2 = p 2 of the plane of the complex variable £ . 
The radius p of this circle is determined in terms of a and h .by 
the formula 



The number p is evidently less than one; h > a. 

The region occupied by the flow is thus transformed into the space 
enclosed between the two concentric circles |i'| = 1 and | £| = p. 

J 

The numbers X and p introduced previously have simple geometri 
cal meanings: X is the length of the tangent drawn from the origin of 

coordinates to the cylinder, and p is the tangent of the quarter 
angle 2 subtended by the cylinder at the origin. 

From formula (l) z is determined: 


z=Xi 


i-K 


(4) 


and from this value are found the values of x and y in terms of the 
argument (n and the modulus r of the complex number c : 
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21r sinu> 

1 — 2rcosu>-^-r , 

v (5) 

y 1 — 2rcos®-f-r* 


‘Skx 

rsm«. = - ^ I+0r + k? 


rcos «> = 


jc , + y s -K* 

* s +( y+W 


( 5 ') 


The transformation of boundary conditions (i) and (il) of the problem 
follow. 


Condition (il) will be considered first. Formula (5) shows that 
this condition may be written as 



1 — 2;>cos»-|-p I 


(no 


In transforming condition (i) it is noted that from formula (5) 
there follows for y = 0 


dy 


Ufa 

1 — COS in 


from which condition (i) assumes the following form: 

[dr 1 — cost»J r = I 


(I') 


where 


=-|>o 


Conditions (II') and (l f ) may be given another form by introducing the 
complex stream function: 


( 6 ) 
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The conditions then assume the following forms : 

C)*^J=o for|C| = I 


(I”) 


Imag.(ar — cz)==« for |C| = p (II") 

The function w is holomorphic at all points of the ring p< |£| < 1, 
hut it will be assumed that the flow around the cylinder has a definite 
circulation. The function w will therefore not be single-valued in 
the ring under consideration. 

In order to obtain the required function w({;), the form of the fol- 
lowing infinite series will be used: 


H“00 

™ = (?) 


with the undetermined coefficients p m and q. 


4. COMPUTATION OF COEFFICIENTS OF SERIES 
REPRESENTING FUNCTION w(£) 

First considered are the relations obtained when condition (II'*) 
is required to be satisfied by the function w(£) . 

In the notation 

Pm = Am J r 

( 8 ) 



the coefficients Am and Buj are real numbers, and T is the circu- 
lation of the flow about the cylinder. 

Condition (II*') then easily transforms into 
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f o»» (H m cos ibw -\- A m sin w«>) — i cR (1 ~! 2, - ; -2,- ; • • • ) * 


m = — ' 

from which the following relations are obtained: 


— In ? 4" £ 0 = LC + a 


-I- -^-2=2 >-f0' 

i /n i r rn 


o'M 


>1-, 


m - m 


m = 1, 2, 3... 


(9) 


( 10 ) 


Condition (I") is now considered. When the following formulas are 
used 

co 

Jmag. (2u«;.] = 2u B n -f- 2 : iJ^ (B m -‘ r B J cos dm +(A m - A_ m )sin mm 

m = 1 

Iniag. + + 

Z < ) 

\ ('» -r i - -2 rnp m -p (/« — l)/>m - 1 ) C m 

m = — :.o 

no difficulty is encountered in representing the condition 




in the form of the fallowing trigonometric series : 


2fvJ9 # + 2p ^ | (B*+B-*)C09;mm ^ (A m — A 

' +OP 

|(m -}-■ 1)&>4 - 1 — 2mBm -}~(/w — — \ j cos -f- 

m — — -,oo 

|(n -j- l)i4»+i.— — 1 )j4*-i] sin »*«> j-==0 


This equation gives a series of relations between the coefficients A m 
and 
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— b 0 — -1 r— 2 |»s, 

B-i) = (B 1 — ,)— (2*,— B- j)— ^ 

+ £_ m ) = m *) ^±J- (£„ + l -B_ ( „ +1) )_ 


m l 


r (-®m - l — fl_ («_ 1)), 


Jm=2, 4. 6 ... } 

^ (/!« - = m (A m 4- ^_ m )_ -^±1- (4„ + ,+* _ ( . + „ )- 


m— 1 
2 


(A m - 1 rf- A _ < ffl _ n ) 


(« = 1, 2. 3. 4, 5, ... ) 


The last of these recurrent relations will he discussed first, 
recurrent relation (10) 


A-m = p tm A m 


as a result of which relation (14) may he written 


^ J 

— 2 — (1 -f P 5 " 1- *) 4»-. + fr(l — P ,m ) — m (1 -f P 5 «)J A m -f 

- 5 ±- L ( | +p ! " +, )^.+. = o ■ 

|«= 1. 2. 3, . . . 1 


For briefness, 


m 

T 


0+P^)4.-V 


• 9 
( 11 ) 
( 12 ) 

(13) 

(14) 
From 
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With these notations the obtained, relation can he rewritten in the fol- 
lowing form: 


*«.-! + S » 1 — 0 

(m- 1,2.3,...! 


(16) 


The coefficient Aq is a certain finite number so that ko = 0 • 

Relation (l3) is now considered. Eliminating from this relation 
the numbers B with negative indices and making use of relation (10) 
yield 

' (1 + P Sm '*) B m -y + 1? (1 - ? tn ) - m (1 4- P 5 ")] + 

- (■ + P’"* 1 ) B.+i = f ] 


For brevity. 






(17) 


is written to obtain for x m the following relation 


m — 1 . m -f- 1 


-1 + S m % m 4* *«.+) SkCp*" | 2pi 


r—m — 1» 


t = 2, 3, 4, 5. . , .) 


If 


(18) 


* 0 — 2« 


(19) 


relation (18) for m = 1 gives condition (12) . 

The system of recurrent relations (16) and (18) is thus obtained 
for the coefficients 


Ai f j4 2 , 

B ly B 2t 


A% , . . . 
B$, . . . 


To these, relations there must he added relations (9) and (ll). 



5. INVESTIGATION OF RECURRENT RELATIONS BETWEEN 


NUMBERS k m 

The relation 

k m-l + s m k n -f k m+1 — 0 

| m = 1, 2. 3, . . .J 

may be used to set up two functional relations between the two 
functions 


«©=£*.?• 

of the complex variable £ . 

First of all 

/ r (C) = 2 V — — 

w m 

or 


( 16 ) 


( 20 ) 


(ff^l 

The new function S(£) is introduced 


OO 



from "which 
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The function S(£) is connected with the function a(£) by the fol- 
lowing obvious relation: 


S(C) + S(p 2 C) = 3 (C) 


(23) 


Together with this relation there is still another obtained from the 
following considerations. Replacing k m on the right side of equa- 
tion (20) by 

— S w-l 


yields 

oo - oo 

m — 2 m — 3 


or 

oo 

(1 + :•-) * C) = *, C - : £ S m _ , ( 24) 

ff*=2 


Transforming the infinite sum on the right side gives 


or 


c b 

5 m — 1 K m— 1 * 


, 1 n-m — 2 r«“l 


5 k rm 1 — 2m* f — P m k C m ~ • — 

Km — ^ I ^7 , 1 -Up-w - "m-l* r 

m= 2 0 


- 1 41_2 o(C) 


The function under the integral may be expressed as follows in terms 
of the functions a(£) and S(£): 



l_.pi*-* 

1 + p*"~* 


k m _ 1 ^-' = 2 5(C)- 6 Q 
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hence 



i'' 


«-»(!) 


Substituting in equation (24) the obtained value of the infinite sum 
produces a new relation between the functions S(£) and cr(£) : 


(i — 9*3(cy— 



25(C)- o(C) 

.c 


d C 


Rewriting this relation in the differential form, 


d r o-y»o 

dr L C 


— — 2j* • 


25(C) -« (;) 
C 


(25) 


and adding to this differential equation the functional equation 


5(C) -{-5(^0 = o(C) 


(.23) 


furnish the two necessary equations for the determination of the func- 
tions cr({;) and S(£). The function S(£) being determined from these 
equations, the function F(£) may be obtained with the aid of rela- 
tion (22) or the relation 


r 



The function 

•CO 

m~l 

may be represented through the function F(£) as follows: 
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M:>=f(4) 

The complex stream function (7) is represented thus: 

w = w,-\- iuK 


Setting 



and 




the function w^(^) may be constructed from the function F(£) 
formula 





or 


w,(0=yi.+ V^(c-+4) 


In the variable z this function may be written 


OO 

M _ A , V A (Z--U)' m + P °"n(z + uy-" 
u j \ z ) — 1 f ^ A m * _j_ ) 3jrn 

m = 1 


The function w-> must now be investigated. 


(26) 


(27) 


(28) 
by the 

. (27-) 
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6. INVESTIGATION OF RECURRENT RELATIONS BETWEEN NUMBERS x m 
In the previously obtained recurrent relation 




, + *«+i = 2>.f P 


2m 


m — 1 



[m= 1, 2, 3, . 




(29) 


the three functions of the complex variable £ may be considered: 


V 


B r . n 


m— 1 


’T0= 


no 

2 


S'G) = 


It 


+ J 


r.2m ’ 


Between these functions there exist the following two easily obtained 
relations : 


5' (C)-fS'(p 2 C) = 3' (C) 


f'q== 2 j s '( :) d ; 

o 

Equation (30) gives the first relation between the functions S'(£) 
and cj'(£).j for obtaining the second relation, the recurrent relations 
(29) are used as follows: • . 

or co 

* G) =* — 2 s «-i *— 2 x * s 4- *i ' + L G) 

m— 2 


(30) 

(31) 


or 


QO 

L (C) = 2U ^ p*« ^ P ! - m - pj C«+> 


(32) 
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Transforming this relation gives 

OO 

(l + ' s ) 3 ' (C) = — C 1 *m 4 + (>)• 


2^s m l x w , c« 1 ;=2 (a j (?) <r, - 2s' o 

f» — 2 O 


(by the computations of the preceding section) . Hence 

r 

( 1 - :y 3 ' (q = - 2/, f 25< ft) d: Xo ,j_ X| ; - f- l q 


from which is obtained 


d ro - :r-o'(C) i 

2yQ- 3 ' ( ') 

. d fMOl 

<* [ •; j “ 

Z J *> • „ 

■ ’ ‘ i f. [ J 

The function L(£) may be represented in finite form: 

1 L Q 

f.*(C — 1)- 0* 

<1 0* ^ 

2 t.c ; 

“ 2(1 — r\y 2 

1 —y\ 

whence 



1 d L C) 

p*(i-p s ) ft-i) 

y ji 

2>f d\ ; 

0 — ?*:)» 

d -p *o s 


Substituting this value of the derivative ~ in equation (33). 

yields ^ 


_d_ (i qfz'c,) 
d’ 


n.. 2S'C)-3'G) r 

■ z . ■ - o- 


o-f^ft-i) ji 

' (I (i 
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This equation together with equation (30) may serve for the determina- 
tion of the functions c r'(0 and S'(0* The function S'(£) being 
found, F’(0 may be computed by formula (3l) . 

The function V 2 of the. preceding section may be written as follows : 


oo 







(35) 


The second sum may be - transformed into 


oo oo 

B_ m C— — 




or 


Thus 


or 


m— 1 


2<*p* 


T— p* 


%P)=s.+ ^-mc+^-+PG)-f'(-f) < 36 > 


* ! (C) = B .+^..nC+ + V *.(«■-£) (37) 


In the variable z this expression may be written 


/ \ o i r t z—u 
w t (z) = B u+-nrln 


2Up{zj±U) 


2 it z + U 1 (1 — p*) 2 T — >-/ (! H-P*) 

n (z — >. if* — p- m (z -f- hi) 9 " 

Zj ‘ (z s -}- 


( 38 ) 
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7 . FORMULA FOR COMPLEX STREAM FUNCTION 

On the basis of formulas (27') and (36), the expression for the- 
function w(£) may now he stated 


® O - (A. 4- * B 0 ) + ~ In ; 4^,- i -\ 

-- - — 9 


FQnFli) 


+ i 


F'Q — F' i K - ') 

\ ’ / j 


(39) 


or 


tc-G) = (A 0 -HB 0 ) 


li 

2" 


lnC4- 


2>.cp- 


{ 


OC CO 

m l m — 1 


, -XV 

m \ r 1 


(39*) 


The number p m is the conjugate complex of the number 

v Pm ~ 

Equations (9) and (ll) permit the determination of the coef- 
ficient Bq and the number a: 


at 





Y — 2|a 


r 

2« 


lnp-f 



1 - 2.1 


l c 


(40) 


The number Aq remains arbitrary, and the circulation T likewise 
remains arbitrary . 

Formula (39*) may be given a more symmetrical form by introducing 
in place of the coefficient p m new coefficients r m and r m by the 
formulas 


r <i=Pm — i >C 

r m =~P m J r ihC 


( 41 ) 
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to obtain 



(39") 


8. APPROXIMATE SOLUTION OF PROBLEM OF DEEPLY SUBMERGED CYLINDER 

The problem of determining the flow of a heavy fluid about a cir- 
cular cylinder is thus reduced to the solution of the two systems of 
the functional equations (23) , (25), (30), and (34). 

Since it is not possible to solve these equations without the aid 
of infinite series or in a finite combination of elementary f un ctions, 
an approximate treatment of the equations of the problem is proposed. 
For this purpose the number p 2 must be considered. In the following 
table are given various values of the ratio a/h and the corresponding 
values of p 2 and p 4 . 


a 

~h 

o 

tr 

*j\ 

0.500 

60" 

0.0718 

0.005151 

.423 

50° 

.0402 

.002416 

-342 

40" 

.031 1 

.000966 

.259 

30" 

.0173 

.000300 

.171 

20° 

.0077 

.000059 

.0*7 

10° 

.0010 

.000004 


This table shows that for a ratio a/h less than 0.342 s l/3, the 

number p 4 does not exceed 0.001, which justifies the rejection of all 
powers of p 2 starting with the second in considering the motion of 
the stream for a/h < l/3. The results here obtained axe in somewhat 
complicated form; therefore only formulas which are suitable for the 
condition at which it is permissible to reject the components with p 2 
are presented. The preceding table shows that this can be done by 
starting, for example, with a/h < l/3. 

By rejecting the terms with p 2 equations (23) and (25) may be 
rewritten as follows: 


d_ 

~cT, 


5 (Cl -- ; (C) 

(1. ') 2 '(C) 
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The function F(£) is defined 



from which 



Equations (30) and (34) are now considered. These, with the preceding 
approximation, may be written as 


d 


j <l-<y3'(Q j 


’'(C) 



r 

2 it 


from which 


V(C) = 


2^ 
1 — 


o-v 




The function F'(^) is now obtained: 



r 1 r- (i-.; ) : 3 ' (C) 
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The complex stream function may now be found, 
mula (39) is used to obtain 


w Q +•«•) + ^ ln : + y 


n 7 

itji ' 


For this purpose for- 


t 1 


+ 


ri 

2™ 


e 


2 ^: c 2*: 



0 


d \ 


(41) 


The constant Bq may be determined by formula (40) . Since = 2x^ 
(by relation (17)), 


2*. 


b ' = T 


_r 

it 

"27" 


In order to determine the remaining constants k]_ and x-j_ the con- 
dition of the absence of disturbance of the flow ahead of the cylinder 
is used. The original variable z in formula (4l) is used: 



H . — n* I 1 r 

~z 4- ki ' 5* 


•* — ki 


n 




In 


z—ki 

z+ki 


*1 4~ V ^j* . 

I* 


-%-i r f i -SLt 
e +LL e * 


w 


dz 


Integrating by parts , 


w 


But 


Q = (.A,+ + B 0 i - — W-^-ln ^41 - A+-V ? *'4. 

\ |» mi J 1 2* z-\-u ji 1 

„ * SL t 

R'-b'C-Z-* 

* J z + h 


z ZL 


- co -f- w 


f *** dz=-A 

4 1 

j z J 

j 


J 




; cos -^4- 2>. sin 

c- C“ 


>./ 


— 00 


4/.* -j-r 


rfS ~ -ie - u I ; tt 

J Z-r hi 


dz 



22 


NACA TM 1335 


Hence 


= {A 0 -f — -f- B 0 i — ) -j- XL in _[_ 

V -«• / 1 2 - z -i- 1 


Vfv ^ 

u 


T 

-fr n f ' C0S 1^ >‘ 2/ sin fi- 
re -f — / - £L 

« J 4/. : -|- $- 


<« 


S' 


2 4- /.? 




The condition at which the function Y becomes zero for y = 0 and 
x = - * will be satisfied if the brackets are equated to zero: 




|tIY 


- 3 \x 




;cos-^-t- 2). sin -^4 
c 2 1 r ! 

4 )*-!-** 


and in addition the imaginary part of the parentheses is equated to 
zero: ' 



The number A q is arbitrary; the real part of the terms in the paren- 
theses is also equated to zero to give 



In this way all the constants have been determined and the complex 
' stream function of the total flow may be written' 


W = — cz — 


2*i 


In 




z -j- >./ 


r — g — 

Ire c * 

TZl 


7 


ef_ 


dz 


(42) 
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Hence for the case here considered of the deep submersion of the cylinder; 
the flow consists principally of the circulation of the velocity about 
the cylinder. The sources of different strengths entering the general 
formula (39") located at the point z = \i begin to be effective only 
at small depths of submersion of the cylinder . 

The vortex representing the cylinder is situated at the point 
z = Xi; this point being somewhat displaced with reference to the 
center of the cylinder. 

Now formula (42) is applied to the computation of the pressure of 
the stream on the cylinder. 

Denoting by & the mass density of the fluid and by X and Y 
the components of the resultant force yields, by the formula of 
Chaplygin, 


Y + iX = 




dW 


The integral is taken over a closed contour surrounding the cylinder: 


dW 

dz 



— OO 


Substituting the preceding in the formula of Chaplygin and applying the 
theory of residues yield 


A' 





c — 



d\' 


CO 


f 


£Z . . 

zcos Sin 



dz 


( 43 ) 


0 



V 
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The first formula determines the wave resistance of the cylinder 
and the second, the lift force. 

For \ = - od the second formula gives the theorem of Joukowski . 

By- formula (4), section 2, it is now possible to find the form of 
the free surface of' the fluid. The equation of the surface of the fluid 
may be written in terms of the function W in the following form: 




t—x+ol 


Making use of formula (42), 


r A cos 


COs4(*— «) — 1 SID 4 (X-S) 


— cb 




— d\ 



Making x approach 


gives 





Hence, far behind the cylinder the surface of the fluid carries steady 
waves of length L = 2nc^/g. 

The amplitude of these waves is 
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For determining the wave resistance, energy considerations are applied. 
The wave resistance R is determined in terms of the energy E of the 
steady wave by the formula 


R 


c — U 

c 


■E 


where U is the group velocity, and 



The energy E is determined by the formula 


£=T e '- A '- 


Making use of all these results formula (43) is again found: 



2gk 

e c '~ 


9. MOTION OF CYLINDER UNDER SURFACE OF 
FLUID OF FINITE DEPTH 

Under consideration is a circular cylinder of radius R moving 
with constant velocity c under the surface of a fluid of depth H. 
The problem is proposed of finding its wave resistance as a function 
of the submersion depth h of the center of the cylinder. This prob- 
lem is solved under the assumptions of Lamb; that is, the magnitude*. 
h/R is assumed to be small and the entire cylinder is replaced' by a 
dipole of moment -ZncBE . 

The motion of the fluid is studied in relation to a system of 
axes of coordinates x and y displaced uniformly together with the 
cylinder. The velocity of the approaching flow at infinity is c. 

By 4>(x,y). is denoted the potential, of the wave velocities. The har- 
monic function $(x,y) satisfies for y = - H the condition 
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Along the free boundary y = 0 the function <i>(x,y) satisfies the fol- 
lowing condition: 


d2<t> 

Hx* 


dx ^ c * dy 


-0 


( 2 ) 


The expression p > 0 is the Rayleigh coefficient of the dissipating 
forces . 

Near the point y = - h, that is, near the center of the cylinder, 
the function $(x,y) on the basis of the above considerations must have' 
the following form: 


<*> = 


cR*x 


* ? + 0 + A ) 5 


Then 


cRT-x 

,,== ”* , +'<y + A)* 


(3) 


and 


<t> 2 = 


, cR*x 
-x' + (y-h? 


(4) 


where $3 is the velocity potential of a certain fictitious dipole 
situated at the point y = h. 

In place of the function <£ a new function <p(x,y) is introduced, 
which is valid for the entire region occupied by the fluid, setting 

$ = <**1 + + ? (5) 


The harmonic function <p(x,y) satisfies the following conditions on 
the boundary of the fluid: 


di _ 

dy ~ dy 


d<t>, 

dy 


for y = — H 


(!') 


<?** ■ „ 1 g d % 2 g d*t 

dx* ' ^ dx 1 c* dy c- dy 


for y = 0 


( 2 ') 



MCA TM 1335 


27 


The function cp(x,y) is obtained in the form of the following definite 
integral: 

00 oc 

cos kx -[- C sin £x)cosh- kydk-\-j ( B cos kx4- -D sin kx) sinh Rydk (6) 

0 


where A, B, C, and D are unknown functions of the variable param- 
eter k. To determine these functions condition (2 ! ) is first considered: 


GO 


] 2 cosh R*x 

~w\y-r wrw 


= eR' 

0 


sin ftxdk 


Substituting this result and (6) in condition (2') yields 


/[( — AM 4- v-kc -f ^ B ) 

e 


cos kx -{- / — k'-C - \JtA -f- 





dk 


OO 


= — R 2 J ke~ hk sin kx dk 


This relation gives the following two equations : 

kA - : iC - - g -B = 0 
‘ r- 

pA 4- kC ArD -~— R* e~ kk 
‘ c- c 


(7) 


Employing condition (l 1 ) gives 


[■ T 2 cR-x(H -A) 

I tv — \*mh f>rr 


cR 2 jke (H h)k sink xdk 

o 


rd<t>- 2cR 2 x{H + h) 


OO 

c/? : J ke~ {H ~ rh)k sin k xdk 

o 
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Applying the results to condition (l') yields 


OO 

f[ 


k (B cosh kH - A sirihkJi) cos kx 4- k(DcoshkH C sinh kH) sin kx dk 


00 

= 2cR- j Ae" ,rt sinh kh sin kxdk 
o 


from which 


B cosh kH - A sinh kH 0. 


D cosh kH — C sinh kH — 2cR-e Hk sinh kh 


or 


B --- A tanh k H 


- D--C tanh kH -L- 2 cR*e~ m - lnh ^ 

' cosh kH 

t 

These equations together with equations (7) permit determining the 
unknown functions A, B, C, and D. 

From equations (l) and (8) the two equations for A and C are 
found: 

* (*_ SL tanh kH)A — \>C = 0 

c 




Solving these equations determines A and C 
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2 g c osh k(H-h) 

c .cosh kH 

l* 8 + (*— 4tanh*/y)* 



cosh/; (// — /t) 
coshkH 


(£—-£- tianh kH) 


\S + (k--^tanh kH? 

C 


Equations (8) then permit finding B and D.' In this way the func- 
tion <p is determined, which makes it possible to obtain the velocity 
potential <t>. The complex potential of the absolute velocities must 
be found: 

W==0-{- V 


First 


0i = Reel { cR ' 


hi 


0, = Reel 


( c * 2 _\ 
V z ~ hi ) 


CO 

9 =• Reel J [(A Di) cos i kz + (C — Bi) sin kz) dk 

0 


where z = x + iy. Hence 


CO 

w ~ h". H- I P -!- Di) cos kz-\-{C — Bi) sin kz J dk 

o 

i 

i 

i 
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The complex potential W of the relative velocities is obtained by 
adding to the right side the term -cz. In terms of this potential 
the pressure of the stream on the cylinder is obtained by the formula 
of Chaplygin: 


X iY 


1 


'dW\- . 

' cj— ’ dz 

\ dz 


where the integral is taken’ over a contour containing the dipole 
z = - ih. 


10 . . COMPUTATION OF WAVE RESISTANCE 

The problem of finding wave resistance is restricted only to the 
computation of X, and therefore only to the imaginary part of the inte 
gral of the formula of Chaplygin: 


* = - 2 '' lm ag- / ) dz 


0) 


The function W may be written 


W; 


cR- 


z-\-hi 


-f G(z) 


The function G(z) is holomorphic about the point z = - hi: 


0(z) =— cz + z C ^:jf+j M + Di) cos kz+(C—Bi) sin kz\ dk 


OO 


( 10 ) 


From formula (9) is obtained 


dW 

dz 


cR 1 

(2-H0 5 


G'(z) 
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from which 



2c/?- 

(z-f-Ai) 1 


G' (z) dz = 4~icR-G”( — ih) 


Therefore 


X = — 2*/?-oc-Reel G" (— ih)) 


Making use of formula (10) the real part of the second derivative of 
the function G(z) at z = - ih is determined: 


CO 

Reel[G"(— ih)] = — j k*(A cosh kh — B sir>hkh)dk 
0 


from which 


X — 2 -R-’jc j k' (A coshhh — B sinh kh)dk 


But, as was found in formula (8), B = A tanh kh. Hence 


X = 2r.R-'jC I k- A C ° s - 

. cosh kH 


Replacing A hy its value gives 


V 1 ' m /■ cosh ’k(H—h) k-dk 

X— 4~',gR*<L j cosE 5 7p>r g— — 

-q p- -f- (R — tanh. Rriy 


which is the formula for ~the wave resistahce in the presence of dis- 
sipative Rayleigh forces . Freed from these forces 
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L lim I n. cosh - k(H-h) k-dk 

: ‘ = °i ' COSh! *" ,= H-(*-^rtanh*W 


Anew variable of integration £ is introduced, setting kH = £ and 
with the following notations: 


h 


\lH »JL J 


In this notation L is rewritten 


CO 



0 


I V 


c osh 1 (1 — a ) £ . £*</* 

’ cosh ’5 V ' ih’+G-xtanhSf 


(11) 


If the number x is less than unity 

f>V£7/ 

then L = G and the wave resistance is equal to zero. When x > 1 
the roots of the equation in E, are considered: 


!V-f-0 — * tanhS) s = 0 ( 12 ) 

This equation may be resolved into the two following equations : 

; = x tanh - -)- t *,t 

S = x tanh ? — <v 

Only the first of these equations will be considered . In view of the 
fact that x > 1, the equation E, =' x tanh £ has_one real root which 
is denoted by E,q. Equation (ll) has one root £, approaching E,q 

as approaches zero; ^ = lim By applying these notations, the 

P=0 

left side of equation (12) is expanded in a series about the point f: 
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!V -I- (: — * tanh t ) 1 (« — •/. tanli* — a,/] (s _ * tanh • -f- a,/] 


= -! A l' d t~ (' — * tanh ; — Ujl). _. s •(; - - ;) -!■ 


(13) 


2 -MM 




COSh 2 ;, 


(5-5)4- 


The integral (ll) may he represented as the sum of a certain contour 
integral taken along the path OIK* and a residue multiplied by 2 ni 
of the function under the integral relative to the point £ = 


I. 


1 

H- 


!im 

S*i = o . 


(OIK) 


-j- lim 

: a i = 0 


2 r.i 

w 


Res | •» 


c osh'O — a) ; 
cosh 2 ; 


(l — x.tanhT ) 2 


The first component on the right side gives zero in the limit. From 
series (13) the limit of the second component can likewise easily be 
found. The following expression is finally found: 


L r • V cosh*(l— a ); 0 
H cosh 2 ;„ — x 


For the expression for X, 


4 

H- 


V cosh 2 (1 — a) ; n 
cosh* ; # — x ' ~ 


(14) 


is now obtained. This is the formula for the computation of the wave 
resistance for velocities c less than the critical velocity -/gH. 
The value £ q is the real root of the equation 


5„ — / tanh ( 15 ) 

From formula (14) is readily obtained the formula of Lamb for the 
wave resistance of a cylinder moving under the surface of a fluid of 
infinite depth. When H is very large, the root £ 0 of equation (15) 
is = x > from which 
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cos h~(l »)Eq , 2j* 

cosh* 5, — * ~ 


Substituting this’ expression in formula (14), the formula of Lamb is 
obtained : 


X 


or 



__ * g,t 

c~ 

e 


Now formula (14) is investigated. In place of £q is written £. In 
the second factor on the right side x is replaced by £/tanh£ to 
obtain 


S'-lgR* tanh s • cosh’-’ (1 — a) l 

H * sinh‘2; — 2: (l^) 


The parameter x varied from 1 to ® so that £ will vary from zero 
to infinity. For £ = 0 the second factor on the right side of this 
formula has a value equal to 3/4. It will be shown that this value is 
a minimum or maximum. Denoting by A the investigated. factor gives 


A — | . ;’sinh 2; cosh{ 1 — a);| 

sinh2; — 25 cosh' j 

Then 

J sin h 2$ 3 ' . , 7 _j_ \ 

2 sinh2$ - 2; 4 ^ 15 5 1 / 


from which 


co sh(l — a ) q 
cosh ; 


=-- 1 a (2 — a);- 






which shows that 
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=0 

l , d\ 



Thus for E = 0 the wave resistance will he an extreme value . It will 
be a minimum if a < 1 - -/8/l5; for larger values of a the wave 
resistance is a maximum. 

For £ = « the factor A becomes zero. This investigation 
explains the general features of the dependence of the wave resistance 
on the velocity. 


11 . CONCLUSIONS 

Formula (16) obtained in section 10 for the wave resistance of a 
cylinder of radius R moving under the surface of a fluid of finite 
depth H led to the following conclusions with regard to the change in 
the wave resistance with velocity of motion and with depth h of its 
submersion: 


The magnitude 



was computed for a series of values of a = h/H: 


* 1 

8 *16* T I 




i.y 

04 ; $ 5 ; 


5 

T? 


a 




8 ^ 


%{ 


for the Froude number f = c / -/gH varying between zero and 1. The 
results of the numerical computations are given in figure 1/ 

The curves of this diagram show a considerable wave resistance 
for small depths of submersion of the cylinder. With increased sub- 
mersion of the cylinder the wave resistance drops sharply for most of 
the values of the Froude number. 

For the value of the parameter a = 1 - -^8/l5, the character of 
the maximum or minimum changes at the point f = 1. For values of 
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a < 1 - -y/8/l5 the curves of the wave resistance have a sharply formed 
maximum. With a increasing from zero to 1 - -Jq/±5 this maximum 
decreases . The presence of this maximum shows up also on the curves 
corresponding to values of a somewhat less than the number 1 - -^8 /l5 
the curves of wave resistance for a between 1 - -y/8/l5 and 0.3 
(approximately) have two peaks, one at the point f = 1 and the other 
near the point f = 0.55. For values of a > 0.3 the wave resistance 
increases monotonically with increasing. velocity j the rate of increase 
of the wave resistance is considerable for Froude numbers f near the 
critical number f = 1. ■ 

In regard to the problem of the first part of the paper (sections 2 
to 8), the results may be. described as follows: When the motion of a 

circular cylinder under the surface of a fluid is studied, the cylinder 
may be replaced by a vortex if the ratio of the radius of the cylinder 
to the depth of submersion of its center is less than l/3 (see table) . 

In the composition of the flow and the lift force the principal part is 
played by the circulation of the velocity about the cylinder. However, 
by considering the general solution of the problem the terms obtained 
by Lamb and Havelock in their work are found. The effect of these 
terms is appreciable only for small depths of submersion of the 
cylinder; for larger depths the circulation of the stream velocity 
plays the fundamental part . 


Translated by S. Reiss 
National Advisory Committee 
for Aeronautics 
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